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Abstract The [O III] A5007 Planetary Nebula Lumi- 
nosity Function (PNLF) is an excellent extragalactic 
standard candle. In theory, the PNLF method should 
not work at all, since the luminosities of the brightest 
planetary nebulae (PNe) should be highly sensitive to 
the age of their host stellar population. Yet the method 
appears robust, as it consistently produces < 10% dis- 
tances to galaxies of all Hubble types, from the earliest 
ellipticals to the latest-type spirals and irregulars. It is 
therefore uniquely suited for cross-checking the results 
of other techniques and finding small offsets between 
the Population I and Population II distance ladders. 
We review the calibration of the method and show that 
the zero points provided by Cepheids and the Tip of the 
Red Giant Branch are in excellent agreement. We then 
compare the results of the PNLF with those from Sur- 
face Brightness Fluctuation measurements, and show 
that, although both techniques agree in a relative sense, 
the latter method yields distances that are ^ 15% larger 
than those from the PNLF. We trace this discrepancy 
back to the calibration galaxies and argue that, due to 
a small systematic error associated with internal red- 
dening, the true distance scale likely falls between the 
extremes of the two methods. We also demonstrate 
how PNLF measurements in the early-type galaxies 
that have hosted Type la supernovae can help calibrate 
the SN la maximum magnitude-rate of decline relation. 
Finally, we discuss how the results from space missions 
such as Kepler and Gaia can help our understanding of 
the PNLF phenomenon and improve our knowledge of 
the physics of local planetary nebulae. 
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1 Introduction 



The [O III] A5007 Planetary Nebulae Luminosity Func- 
tion (PNLF) has been a reliable and precise ex- 
tragalac t ic dis t ance indicator f or ove r ^ 20 years 
( Jacobvl Il989t ICiardullo et al.l Il989at Ijacobv et al. 
19921 ). During this tir ne, the method ha s been ap 



plied to both ell i ptical (jJacobv et al.lll990l) and spiral 
Feldmeier et al.l 19971) galaxies, to galactic spheroids 



Ciardullo et al 



11989a: 



Jacobv et al 



Hui et al. 



- -- - - 19891 : 

19931 ) and disks (JFeldmeier et al.l Il997l ). and even to 



stars lost within the intergalac tic environment of rich 
clusters JFeldmeier et al.l I2OO4I ). In addition, PNLF 
observations are relatively easy: the method requires 
neither space-based observations nor heroically long in- 
tegrations, and the photometric procedures needed to 
derive accurate distances are simple and straightfor- 
ward. As a result, the PNLF is an integral part of the 
extragalactic distance ladder, and perhaps the best tool 
we have for examining systematic differences between 
Population I and Population II distance methods (see 
Figure [T|) . 

Of course, no technique is perfect. The PNLF is not 
well-suited for small galaxies which contain few plane- 
tary nebulae, nor can it reach to ~ 100 Mpc, where the 
unperturbed Hubble Flow dominates. For all practical 
purposes, PNLF observa t ions a re li mited to ~ 20 Mpc 



(but s ee iGerhard et al.l (|2005l ) and IVentimiglia et al.l 



( 20111 ) for PN measurements in Coma and Hydra!). 
The technique also requires the use of a narrow-band 
[O III] A5007 fiher (FWHM - 50 A), whose properties 
in the converging beam of a telescope are well known. 
In this age of fast, extremely wide-field imagers, this 
can be a severe limitation, especially for systems where 
the 5007 A line is redshifted out of the bandpass of a 
rest-frame [O III] filter. Finally, the PNL F cannot easily 



be cal ibrat ed in the Milky Way (bu t see iMendez et al 



( 19931 ) and iKovacevic et al.l (|2011b[ ) for attempts at a 
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Fig. 1 The extragalactic distance ladder. The blue boxes 
show techniques useful in star-forming galaxies, the red 
boxes give methods best suited for Pop II systems, and the 
brown boxes represent geometric measurements. Since the 
PNLF is equally effective in all stellar populations, it is ideal 
for cross-checking the results of different methods. 



local calibration), and wc arc far away from having a 
theoretical understanding of the method. 



2 Why the PNLF Can't Work 

The universe constantly surprises us. For example, in 
the 1980's, if one were to ask whether 4-m telescopes 
could detect [O III] A5007 emission from planetary neb- 
ulae in the Virgo Cluster core, the answer would have 
certainly been no. The main-sequence turnoff of an 
elliptical galaxy is ~ ^Mq, which, through the initial 
mass-final mass relation, means that the central stars 
of the PNe currently being p roduced have relat ively 
small masses, i.e., ~ 0.52Mq (JKalirai et al.l 120081 ). It 



is well established that cores in thi s mass range cannot 
be more luminous than ^^ 1000 L0 (|Schoenberneiill983 : 
Vassihadis fc Woodlll994l) . and that less than ^ 12% of 



this flux c an be reprocessed into the [O III] A5007 emis - 



sion line ( Dopita et al.lll992t ISchonberner et al.ll2010l ). 



Thus, at maximum, the PNe inside the elliptical galax- 
ies of Virgo should have [O III] monochromatic fluxes of 
1.2 X 10^^'' ergs cm^^ s^^ (^5007 ^ 28.5), well below 
the threshold for night-long observations with a 4-m 
class telescope. 

Of course, one might argue that real galaxies contain 
a mix of stellar populations and that there will always 
be some PNe produced by higher mass stars. But this 



cl aim poses another problem. As elegantly pointed out 
bv lRenzini fc Buzzoni I (|l986[ ). the number of PNe pro- 



duced by any single-aged stellar population is linearly 
proportional to the luminosity of that population, i.e.. 



N{PN) ^ B ■ t ■ L 



(1) 



where t ^ 500 yr is the lifetime of a PN in its [O III] 
brigh t phase (JMarigo et al.l l2004l : ISchonberner et al. 



2007|) , L is the total integrated (bolometric) luminosity 
of the population in question, and B is the population's 
luminosity specific stellar evolutionary flux. Remark- 
ably, B does not depend on age, metallicity, or initial 
mass function: to within '^ 10%, all populations older 
than ~ 1 Gyr have B ~ 2 x 10"^^ stars yr~^ L^^ 



(jRenzini fc Buzzoni 1119861 ). Consequently, in order for 
a Virgo elliptical to have a sizeable number of [O IIl]- 
bright PNe, ^ 10% of its stars must be no older than 
~ 1 Gyr. Such a result has long been ruled out via in- 
tegrated light spe ctroscopy (e.g.. IVazdekis et al.lll997 : 



Trager et al.ll200Gl ) 



Furthermore, even if the PNe were detected, there 
is absolutely no reason to expect the bright end of the 
PNLF to be a standard candle. From stellar evolution 
theory, the [O III] A5007 magnitude of the brightest PNe 
produced by a simple (single age) stellar population 
should change with time, following 



dM; 



5007 



dt 



dh'h 



007 



(ilogL5007 

dlogL* 



d log ^5007 

d log L, 
f dnif 



drat 
diTif J \ dmt J \ dt 

(-2.5).(^l).(^7)-(0.11).(0.8i-i-^) 



-1.5t" 



(2) 



where L* is the luminosity of the central star, m/ is 
the mass of the PN core, mt is the turnoff mass of the 
stellar population, and t is the age of the population 
(measured in Gyr). All these derivatives are reason- 
ably well known via our knowledge of main sequence 
and post-AGB stellar evo l ution (e.g.. Ilben fc Laughlinl 



M 



mass relation 



le.g., 



1989t IVassiliadis fc Wood 1994 ), the initial mass-final 



Kalirai et all 120081: iDobbie et al 



20091) ■ and nebular p hysics (e.g.. iFerland et al.lll998 : 
Perinotto et al.l 120041 ). (The most uncertain term in- 
volves the reprocessing of stellar luminosity into line 
emission at 5007 A, bu t models such as those by 
Schonberner et al.l ( 20101 ) suggest that this derivative 
is a slowly varying quantity.) For old stellar systems 
(t -- 10 Gyr), equation Q implies that the PNLF cutoff 
should fade by ^ 0.1 mag per Gyr, and over a timescale 
of -^ 10 Gyr, th e decline in brightness should be more 
than ^ 4 mag (JMarigo et al.l l2004t ISchonberner et al. 



20071) . In this age of precision eosmology, such a strong 
dependence would have been certainly been seen. 

Finally, an examination of planetary nebula mor- 
phologies in the Milky Way reveals that most bright 
PNe are asymmetric, and have position-dependent dust 
opacities within their circumstellar envelopes (e.g.. 



Ueta et al.ll2000l: ISiodmiak et al.ll2008D . The effect this 



dust has on the emergent [O III] A5007 emission is sig- 
nificant: although each PN is different, the average 
[O IIl]-bright planetary n ebula is self-extincted by more 
than 0.6 mag at 500 7 A (JHerrmann fc Ciarduliol 12009 : 
iReid fc Parkerll2010l ). In fact, as Figure [2] shows, the 
bright-end cutoff of the PNLF is largely defined by the 
effects of sight-line dependent circumstellar extinction. 
At face value, this fact alone would seem to preclude 
the PNLF from being an effective distance indicator. 
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Fig. 2 The [O III] A5007 to Hq-|-[N II] line ratio plotted 
against [O III] absol ute maRnitude for a complete sample 
of PNe in the LMC (JReid fc Parkeij|2010l l. The light blue 
points show the intrinsic line strengths; the red points il- 
lustrate the line intensities that are observed. The dotted 
green line displays the observed value of the PNLF cutoff, 
Af*. The plot suggests that A'l* is largely defined by the 
action of circumstellar dust. 



3 But the PNLF Does Work 

Despite these arguments, the PNLF is a precise ex- 
tragalactic standard candle, capable of generating dis- 
tances to better than ^ 10% in clusters as far away as 
Virgo and Fornax. To demonstrate this, one need only 
look at our Local Group neighbor, M31. As illustrated 
in Figure El the PNLF of M31's bulge {R < 1 kpc) 
has the same bright-end cutoff (defined via the abso- 
lute magnitude M*) as that of the system's inner disk 
(6 < i? < 10 kpc) and outer disk/halo {R > 15 kpc). 
Given the range of stellar population being observed, 
this invariance is remarkable. 



B 

Q. 







20.5 21 21 .5 

Apparent >.5007 Magnitude 

Fig. 3 The observed planetary nebula luminosity func- 
tions for samples of M31 PNe projected at three different 
galactocentric radii. The curves show the best-fitting em- 
pirical law. Although the stellar populations of M31's inner 
bulge (R < 1 kpc), inner disk (6 kpc < i? < 10 kpc), and 
outer disk/halo (i? > 15 kpc) are vastly different, the de- 
rived PNLF distances to each region are consisted to within 
~ 0.05 mag. 



M31 is not the only place where the PNLF has 
been tested. PN studi es in galaxies as diverse as MB S 



( Giardulloet all 2004), NGC 5128 ( Hui et al.Ml993 ) 



and NGC 4494 (jjacobv et al.l 119961) . have been unable 
to detect any change in M* with galactocentric ra- 
dius. Similarly, observations of galaxies in groups an d 



clusters, such as Triangulum (ICiarduUo et al. 



1997; CiarduUo et al.l 



1991 L 



the M 81 Gro up Ijacobv et al.l Il989l : ICiardull o et al 
20021) ■ Leo I dCiarduUo et al.l Il989b[ LFeldmcicr et al 



20021) ■ Virgo djacobv et al.lll990l: 
Feldmeier et al.l l2007l ). and Fornax (JMcMillan et al.l 



1993HTeodorescu et al.ll2005l:lFeldmeier et al.ll2007l) fab 
most) always place all the galaxies comfortably within 
the typical group diameter of ~ 1 Mpc. (The lone ex- 
ception is in Virgo, where the PNLF clearly resolves the 
M84/M86 system which is known to be infalling into 
Virgo from behind.) The multitude of internal tests 
place strong constraints on the types and amplitudes of 
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Fig. 4 PNLF distance measurements in tlie Triangulum 
Spur, the Leo I Group, tlie Virgo Cluster, and the Fornax 
Cluster. The M84/M86 system is also marked: although 
projected onto the Virgo Cluster core, this group is actually 
falling into Virgo from behin d. The Hubble types are from 
ISandage fc TammannI (|l98ll ) . The diversity of Hubble types 
illustrates the robustness of the PNLF method. 



any systematic errors that might be associated with the 
technique. The tests also demonstrate that the PNLF is 
capable of generating relative extragalactic distances to 
just a few percent in a variety of galactic environments. 



4 The Calibration of the PNLF 

There is no theoretical explanation for the constancy 
of the PNLF, and reliable distances to Galactic plan- 
etary nebulae are few and far between ( Phillipa 12004 : 
Harris et al.l 120071) . Hence, the only way to obtain an 
absolute measurement of M* and look for external er- 
rors in the PNLF method is to survey extragalactic sys- 
tems that have distances known from other methods. 



CiarduUo et al.l (J2002I ) did this via PN observations in- 
side 13 galaxies with Cepheid distances derived by the 



(Il989a[ ) 



HST Key Project (IFreedman et al.ll200ll). By adopting 
the "universal" PNLF first proposed bv lCiarduUo et al.l 



} 



(3) 



the authors obtained a value of M* 



-4.47, where 



M, 



5007 



= -2.51ogF5oo7- 13.74 



(4) 



and i^5oo7 is the [O III] A5007 flux in ergs cm~^ s~^. 

Remarkably, despite ~ 20 years of observations, 
non-parametric statistical tests such as Kolmogorov- 
Smirnov and Anderson-Darling still have not found 
any reason to reject the shape defined by this simple 
analytic function (bu t see the slitlcss spectroscopy of 
Sambhus et al.l ( 20061 1 for a curious result in the el- 
liptical galaxy NGC 4697). Of course, this does not 
mean that other formulations of the PNLF are not 
possible. For exampl e, the recent distanc e determina- 
tions to NGC 4697 JM cndcz et aLllJOol. NGC 134 4 



( Teodorescu et al.l [2005 ). M8 2 djqhnson et al.l [2OO9I) 



NGC 821 dTeodorescu et all l2010l) . and NGC 4649 
(iTeodorescu et al.ll2011 ) all use a numerical form of the 
PNLF derived by combining models of post-AGB evo- 
lution with empirical constraints on the excitation of 
[O IIIJ A5007 in Galactic and M agellanic Cloud plane- 
taries (JMendez fc Soffneiill997| ). These measurements 
are self-consistent within themselves, and usually pro- 
duce distances that are very similar to those found from 
equation ([3]). However, because the functions are differ- 
ent, there may be small (< 0.1 mag) systematic offsets 
between the distances derived in this way and those in- 
ferred from the purely empirical PNLF. The quoted 
uncertainties of the individual distances may also have 
slightly different meanings. 

Since 2002, a few of the PNLF and Cepheid mea- 
surements have improved, and the number of galaxies 
studied by both methods has increased to 16. But qual- 
itatively, little has changed. As the left-hand panel of 
Figure [5] illustrates, all stellar populations more metal- 
rich than the LMC appear to have the same value for 
the PNLF cutoff, M* ^ -4.5. In the smaller, metal- 
poor systems, M* does fade, in ge neral agreeme n t with 
the nebular model pre d iction s of iDopita et al.l (|l992[ ) 



and Schonberner et al 



( 20101 ). However, since these 
low-mass systems have few PNe and poorly defined PN 
luminosity functions, this systematic shift in M* is not 
an important limitation of the method. 

We do note that the best-fit value of M* can be 
changed slightly, depending on how one models the 
response of the Cepheid period-luminosity relation to 
metallicity. The solid points displayed in the left- 
hand panel of Figure [5] adopt the Cepheid distances 
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Fig. 5 Measurements of M* derived using distance estimat es from Cepheids (left panel) and the Tip of the Red Giant 
Branch (right panel). Most of the Cepheid distances come from lFreedman et alj (|200ll ) and assume no metallicity correction. 
The open circles show the same calibration with an assumed value of 7 = —0.2 for the dependence of the Cepheid period- 
luminosity relation on metallicity. Such a correctio n increases M* b y ~ 0. 07 mag, whil e incre asing the dispersi o n onl y 
slightly. The TRGB distances come principally from lDalcanton et all (|2009l ). [jacobs et al.l (|2009l ). and lHislop et all (|201ll ). 



of iFreedman et al.l (|2001l ) and assume that the period- 
luminosity relation for Cepheids has no dependence on 
metallicity. If we exclude the small, low-metallicity 
galaxies (i.e., objects with 12 -I- log 0/H < 8.45), these 
data imply a most likely value for the PNLF cut- 
off of M* = -4.46 ± 0.05 (standard deviation of the 
mean), and an external scatter (a = 0.16 mag) that 
is fully consistent with the inte rnal errors of the mea- 
suremcnts. Alternatively, if the lFreedman et al.l ( 2001 ) 
period-luminosity data are used with their suggested 
metallicity dependence (7 = —0.2), then the most likely 
value of M* brightens to M* = -4.53 ± 0.04, while the 
external scatter of the individual M* measurements in- 
creases only marginally to 0.18 mag. Finally, an ex- 
treme value of M* = —4.68 ± 0.09 can be obtained 



if the Cepheid distances and mctallicities of ISaha et al 



(20061) are used in the analysis. Under this assumption 
however, the internal errors of the individual PNLF 
measurements are no longer consistent with the ob- 
served galaxy-to-galaxy scatter, which explodes to al- 
most 0.3 mag. 

Rather than relying solely on the Cepheid period- 
luminosity relation, we can take another route to the 
PNLF zero point by using distances inferred from mea- 
surements of the Tip of the Red Giant Branch (TRGB). 
Thanks mostly to HST, 18 PNLF galaxies of vary- 
ing Hubble types now have reliabl e distances from this 
technique (IDalcanton et al.l l2009t Ijacobs et al.l 12009 : 
Hislop et al.ll201ll ). The values of M* based on these 



measurements are shown in the right-hand panel of Fig- 
ure m If we exclude the results from the two distant 
systems (NGC 3368 and 3627) where the TRGB and 
Cepheid distances are in conflict, the plot looks remark- 
ably similar to that obtained from the Cepheids. As 
before, there is evidence for a decrease in M* in low- 
luminosity, low-metallicity systems. However, in sys- 
tems more metal-rich than the LMC, M* again has a 
best-fit value of M* — —4.53 ± 0.06 and there is no ev- 
idence for a metallicity dependence. The fact that this 
number, with its small dispersion, agrees with the cal- 
ibration derived from Cepheids supports the argument 
that the zero point of the PNLF system is secure at the 
~ 5% level. 



5 The PNLF and Surface Brightness 
Fluctuations 

Both the Cepheid and TRGB methods establish the 
zero point of the PNLF distance scale to ~ 5%. 
Wc can now use this fact to investigate the calibra- 
tion of other methods, and search for systematic er- 
rors in the distance ladder. For example, the Surface 
Brightness Fluctuation (SBF) method is a popular way 
of obtaining high-precision distances to galaxies with 
smooth l uminosity profiles, i.e., elliptical and lenticula r 
systems JTonrv fc Schneider] Il988t iTonrv et al.ll200ll) 
The t echnique is efficient, well calibrated (JTonrv et al 
200 ll iMei et al.l l2005l : iBlakeslee erahl l2010t ). and has 
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Fig. 6 A histogram of tlie differences between tlie PNLF 
and SBF distance moduli for 33 galaxies measured by both 
techniques. The two worst outliers are the edge-on galaxies 
NGC 4565 (A/i = -i.l4) and NGC 891 (A/i = +0.56). 
NGC 4278 is also an outlier (A^ = -0.97). The curve 
represents the expected dispersion of the data. The figure 
demonstrates that, although the internal and external errors 
of the methods agree, the absolute scales defined the two 
techniques are in conflict. 



a solid theoretical foundation (e.g.. IlIu et a l. 2000"; 



Blakesleeet al.ll200lt iMarin-Franch fc Aparicio. .2006.1 . 
The most direct way of obtaining the zero point of this 
method is to use the seven Cepheid galaxies with bulge 
populations large enough to be measured with the SBF 
technique. If the zero point derived from these galaxies 
is accurate, then a comparison of PNLF and SBF dis- 
tances should show good agreement, with a mean near 
zero, and a scatter that is representative of the internal 
errors of the methods. 

Figure IHl displays the difference in distance moduli 
for the set of 33 PNLF galaxies with SBF /-band m ea- 
surements ( Tonrv et al.ll2001t iBlakeslee et al.ll2010l) . It 
is immediately obvious that the systematic offset be- 
tween the two methods is not zero, nor is it close 
to zero. The SBF distance scale is, in the median, 
A/x = 0.33 larger than PNLF distance scale. More- 
over, this property is not restricted to the groundbased 
/-band dataset. Twelve PNLF galaxies have high- 
quality SBF measurements in t he z-band fro m obser- 



vations with the ACS on HST (JBlakeslee et al.l 120091 ): 
the offset between their PNLF and SBF distance mod- 
uli is A/i = 0.42. Similarly, 16 galaxies have PNLF 
and SBF //-band data from HSTs NICMOS camera 
The median offset between these 
0.29. Clearly, there is a problem 



(| Jensen et al.l 120031 ). 

two datascts is A/i = 
with the zero point of one (or both) of the techniques. 
The fact that the offset is due to a zero point er- 
ror, rather than a systematic trend with stellar pop- 
ulation or distance can be shown in two ways. First, 



the solid curve in Figure [6] is not a fit to the data: it 
is instead, the expected scatter in the measurements, 
as determined by propagating the uncertainties asso- 
ciated with the individual PNLF and SBF distances 
and Galactic reddening. The agreement between the 
curve and data proves that the quoted uncertainties 
of both methods are reasonable, and that there is lit- 
tle room for additional random errors associated with 
the measurements. Similarly, if wc plot the individual 
differences in distance moduli against against external 
parameters such as galactic absolute magnitude, color, 
distance modulus, or number of bright PNe (within 
0.5 mag of M*), we find only one significant trend. As 
Figure [7] illustrates, the offset between the PNLF and 
SBF correlates slightly with SBF distance modulus, in 
the sense that distant galaxies have smaller PNLF dis- 
tances than expected f rom their fluctuation magnitude. 
Ferrarese et al.l (|2000l ) interpreted this trend as being 



due to a systematic error in the PNLF that only ef- 
fects galaxies with distances greater than 10 Mpc. At 
the time, such an error could plausibly have been due 
to the existence of foreground intracluster stars in the 
clusters of Virgo and Fornax. However, in the past few 
years, this explanation has become untenable. Not only 
is the hypothesis at odds with the informat ion provided 



by deep Vir go Cluster surface photometry (JMihos et al. 



20051 120091 ) ■ but additional data has shown that the 



distance discrepancy is not restricted to galaxies in the 
centers of clusters. Objects such as the isolated ellipti- 
cal NGC 821 and the x-ray faint elliptical NGC 4697 (in 
the Virgo southern extension) , display the same curious 
offset. 

In fact, a closer look at Figure [7] shows that the ap- 
parent correlation with SBF distance modulus is not 
driven so much by the measurements of the most dis- 
tant systems, as it is by the behavior of the relatively 
nearby spiral galaxies. If these later-type objects are 
excluded from the analysis, the correlation with dis- 
tance disappears. Moreover, except for the systematic 
offset between the two techniques, the only consistent 
pattern visible in Figure [7] involves the dispersion: the 
scatter between the PNLF and SBF distance estimates 
for spiral galaxies is much larger than what is seen for 
elliptical and lenticular systems. 

So where does the offset come from? As Figure [5] 
demonstrates, the PNLF zero point is well calibrated 
by both Cepheid and TRGB measurements; while a 
small increase in the distance scale may be possible, any 
increment greater than ^0.1 mag is excluded by the 
data. Likewise, an examination of the SB F zero point 



calib r ation shows little room for error (e.g.. lAihar et al. 



200 It iFerrarese et al.l 120001 ). The median offset for the 



seven galaxies with both SBF and Cepheid distance de- 
terminations is A/i = —0.10 if the period-luminosity 
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Fig. 7 The diflFerence between SBF and PNLF distance moduli plotted against galactic absolute magnitude, distance, 
color, and number of PNe within 0.5 mag of M* . The red points represent elliptical and lenticular galaxies, the orange 
points are gas-rich lenticulars and/or recent mergers, and the blue points are normal spiral galaxies with SBF measurements 
in their bulge. The correlation with SBF distance modulus, though significant (P ^ 0.01), disappears if spiral galaxies are 
excluded from the analysis. Note that the scatter for spiral galaxies is much larger than that for the early-type systems. 



relation is independent of metal abundance, or -1-0.08 if 
7 = —0.2. It should be noted, however, that the uncer- 
tainties associated with several of the individual fluctu- 
ation magnitudes used in this calibration are large, and 
a similar analysis with the 12 SBF galaxies with TRGB 
measurements yields an offset of A/i = —0.17. Thus, 
while there is some evidence to support a modification 
in the SBF zero point, any downward adjustment to 
its distance scale must be slight, < 0.15 mag. Another 
explanation is therefore needed to reconcile the PNLF- 
SBF discrepancy. 

The answer to this paradox is internal extinction. To 
calibrate an extragalactic standard candle, one needs to 
measure the apparent brightness of the candle, to, and 
adopt values for the distance modulus to the galaxy, 
/x, and the intervening extinction, E{B — V). In other 
words, 



Here, M is the derived absolute magnitude of the stan- 
dard candle in question and R\ is the ratio of total 
to differential extinction at the wavelength of interest. 
(For [O III] A500 7 observations of plane tary nebulae, the 
extinction law of ICardelli et al.l (jl989l ) gives R\ =3.5.) 
For most methods, including the PNLF, if the redden- 
ing to a galaxy is underestimated, then the brightness 
of the standard candle will be underestimated, and the 
distance scale implied by the observations will be under- 
estimated. However, in the case of the SBF technique, 
there is a non-negligible color term, which causes redder 
stellar populations to have dimmer fluctuation magni- 
tudes. Consequently, the absolute fluctuation magni- 
tude, M, is given by 



M, 



C + a(l/-/)o 



(6) 



M 



^l-RxE{B~V) 



(5) 



where a is the slope of the color dependence. The zero- 
point of the SBF distance scale, represented by the con- 



stant C, is therefore defined tlirougli 

C^fni^fi-a{V-I)ohs + {aRv-ia + l)Ri)E{B-V) 

(7) 

For the /-band SBF measurements, a — 4.5, Ry — 3.04. 
and Ri = 1.88, so 



C ^mi -fi + 3ME{B - V) 



(8) 



(For the HST z-band photometry, the coefficient for the 
reddening term is +1.71, while for the IR fluctuation 
data, the value is +5.35.) So, while the response of 
SBF measurements to reddening has roughly the same 
amplitude as that of the PNLF, it has the opposite 
sign! Even if both analyses were performed perfectly, 
our imperfect knowledge of foreground extinction would 
cause the two distance methods to produce different 
results. 

Of course, simple random uncertainties in the esti- 
mates of foreground reddening, such as those associ- 
ated with the extinction maps of ISchlegel et al.l (|l998l ). 
should not produce a systematic error between the two 
distance scales. However, the calibrators for the SBF 
and PNLF techniques are mostly mid-type spiral galax- 
ies, while the vast majority of the methods' program ob- 
jects are elliptical and lenticular systems. This popula- 
tion difference can create a systematic bias between the 
two systems. If the calibrator galaxies have, on average, 
as little as E{B — V) ^ 0.02 more internal reddening 
than the early-type systems that are the methods' main 
targets, the result will be a 0.15 mag offset between the 
PNLF and SBF distances, in exactly the direction that 
is observed. This seems quite plausible, especially since 
the bulges of most mid-type spirals exhibit more evi- 
dence for dust than do the inter iors of normal elliptical 
galaxies (jWindhorst et al.ll2002l ). In fact, given the ex- 
tremely strong amplitude of the extinction dependence, 
it would be surprising if there weren't an offset! 

From the above arguments, it seems clear that the 
offset between the PNLF and SBF distance scales 
is likely the result of several small effects: a slight 
(< 0.1 mag) underestimate in the brightness of the 
PNLF cutoff, a slight (< 0.2 mag) overestimate in the 
zero point of the SBF method, and a small shift in 
both the PNLF and SBF caused by the systematically 
larger dust content of spiral bulges. Without the cross- 
checking ability of the PNLF, such subtle effects could 
easily go unnoticed. 



supernovac in the past century, and about half of these 
objects have well-observed light curves. Included in this 
sample of well-studied supernovae are 6 SN II-P, 2 SN 
ILL, 1 SN lib, 1 SN Iln, 1 SN lb, 1 SN Ic, and 13 SN 
la (sec Table [T]) . The latter can be used to derive an 
estimate of Hq. 

Table 1 SN la m PNLF Galaxies 



SN la Host 



Comment 



1972E 


NGC 5253 


Photoelectric Photometry 


1980N 


NGC 1316 


Photoelectric Photometry 


1981D 


NGC 1316 


Photoelectric Photometry 


1986G 


NGC 5128 


High Internal Extinction 


1989B 


NGC 3627 




1991bg 


NGC 4374 


Subluminous 


1992A 


NGC 1380 




1994D 


NGC 4526 




1998bu 


NGC 3368 




2006dd 


NGC 1316 




2006mr 


NGC 1316 


Subluminous 


2007on 


NGC 1404 




2011fe 


NGC 5457 





Figure [8] displays the 1^-band maximum magnitude- 
rate of decline relation for a set of well-observed su- 
pernovae in the Hu bble flow, i.e., in th e redshift range 
0.01 < z < 0.08 iHickenet al.ll2009l) . under the as- 



sumption that Hq = 70 km s~^ Mpc~^. Also plot- 
ted are the maximum absolute magnitudes of 13 well- 
observed local SN la, using distances derived from the 
PNLF. The agreement between the two sets of data 
is surprising good, especially when one considers that 
several of the supernovae are from the era prior to CCD 
measurements and another is affected by high internal 
extinction. This consistency demonstrates the viability 
of the PNLF as a SN la calibrator, and its utility. Over 
half the local supernovae plotted in the figure occurred 
in early-type galaxies with no significant Population I 
component. Cepheids cannot be used to measure the 
distances to these objects. In order to study the bright- 
nesses of SN la across galaxy types, one needs a method 
such as the PNLF, as it is the most direct way to link 
measurements within Population I and Population II 
stellar systems. 



7 Prospects for the Future 



6 The PNLF and SN la 

Over 50 galaxies currently have robust PNLF distance 
measurements. These galaxies have hosted ^ 50 known 



The greatest impediment to improving the PNLF dis- 
tance scale is the lack of an adequate theory to explain 
why the technique is so robust. The constancy of the 
PNLF cutoff cannot be explained by single-star stellar 
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Fig. 8 The absolute V magnitudes of 13 SN la located 
within galaxies with PNLF distances, plotted against the 
amount of fading which occurred in the first 15 days after 
maximum. For comparison, a sample of supernovae in the 
redshift range 0.01 < z < 0.08 is also shown, under the as- 
sumption of Ho = 70 km s"^ Mpc"^ Except for SN 1972E 
(in NGC 5253, shown as an open circle), all the supernovae 
are from the past 30 years. 



evolution: the energy emitted by an M* planetary neb- 
ula dem ands that the PN has a co re mass of at least ^ 
0.6Mq IVassihadis fc Woodl[l993) . and old stellar pop- 



ulations simply cannot produce such high-mass objects 
The progenitors of '~ O.6M0 white dwarfs are main- 
seque nce stars with masses M 



2M(7, (Kalirai et al. 



20081) . and stars such as thes e have very short mai n 
sequence lifetimes (r < 1 Gvr: llben fc Laughlinlll989f ). 



While early-type galaxies m ay contain a small num- 
ber of intermediate-age stars (jKavirai et al.ll2007l ). this 
population is not nearly large enough to produce the 
number of PNe actually observed. Consequently, in or- 
der to explain the PNLF of elliptical galaxies, a more 
exotic channel of stellar evolution must be considered. 
Binary evolution is the likely solution to this prob- 
lem. There is now considerable evidence to suggest 
that a signifi cant fraction of PNe form via binary star 
inter actions dSokeilbood iMoe fc De Marcol [2OO6I) . In- 



deed, ICiardullo et al.l (|2005[ ) has proposed that an im- 
portant channel for the creation of an [O IIl]-bright 
PN involves co nservative ma ss transfer whilst on the 
main sequence (JMcCrealll964[ ). and the subsequent cre- 
ation of a blue straggler star. Though this hypothesis 
is difficult to test, clues as to its viability may soon 
be pro vided by astrom e tric an d photometric space mis- 
sions. IClarkson et al.l (|2011[ ) have shown that space- 



based astrometry of the Galactic bulge can provide a 
census of blue stragglers in a coUisionless ste llar envi- 
ronment with a well-studied PN population ( PottaschI 
1990t iKovacevic et al.ll2011al ). By comparing the PN 



and blue stra ggler number densities to their expected 
lifetim es (e.g.. lSchonberner et al.ll2010tlLombardi et al. 



20021) ■ one can test whether the stellar merger hypoth- 
esis is self-consistent. Similarly, if blue straggler evolu- 
tion is an important channel for populating the bright 
end of the planetary nebula luminosity function, it must 
also be a significant contributor to the A-star popu- 
lation of the solar neighborhood. Asteroscismological 
studies, s uch as those being performed with the Kepler 
satellite (JBalona et al.ll2011l ). may be able to disentan- 
gle the two evolutionary scenarios, thereby shedding 
light on the creation rate of coalesced objects. 

Another advance that we can look forward to is an 
improved distance scale to PNe in the Milky Way. Most 
Galactic planetary nebulae have distance estimates 
based on such dubious methods as the amount of fore- 



ground extinction (JGiammanco et al.ll2011l ) or the sta- 
tistical correl ation between nebular emission and nebu- 
lar si ze (e.g.. IShklovskvlll956t Ivan de Steene fc Zijlstral 
19951 ). Factor of ~ 2 errors in these estimates are 
not uncommon, even for bright, well-observed objects, 
and as a result, obtaining a loca l calibration of the 
PNLF cutoff is extreme ly difficult JMendez et al.lll993 : 
Kovacevic et al.ll2011b[) . This will soon change, how- 
ever, as Gaia will expand the number of PNe with re- 
liable parallax measurements from ~ 20 (IHarris et al. 
20071 ) to more than ^ 200. This should facilitate dozens 



of studies of the local PN population, including its 
[O III] luminosity function. 

Unfortunately, a Gaia-based luminosity function, by 
itself, will probably not enable us to significantly im- 
prove the calibration of the PNLF. As Figure [5] il- 
lustrates, M* is highly dependent on the amount of 
extinction produced by circumstellar dust, and, de- 
spite the extraordinary progress made in undcrstand- 
ing the physics of the p lanetary nebula phenomenon 
( Schonberner et al.ll2010l ). this component of the prob- 
lem has yet to be modeled. This is major issue, since 
observations of Galactic PNe are affected not only by 
the actions of circumstellar dust, but by foreground 
extinction as well. The problem of disentangling the 
two reddening components may ultimately limit any 
attempt to calibrate the planetary nebula luminosity 
function in the Milky Way. 



8 Conclusion 

The PNLF is neither a primary distance indicator nor 
a method that can reach out of the Local Superclus- 
ter into the Hubble Flow. It has no theoretical basis, 
and next generation missions such as Gaia will not sub- 
stantially improve its calibration. Yet the technique re- 
mains an important part of the extragalactic distance 
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ladder: it is, perhaps, the best tool we have for link- 
ing the Population I and Population II distance scales, 
and identifying systematic offsets between the different 
methods. Moreover, it is still a relatively efficient tech- 
nique that does not require space-based observations or 
extremely long exposure times. Since PNLF distances 
are ofte n produced as a by-product of kinematic stud - 
ies (i.e.. [Herrmann et al.ll2008t iTeodorescu et al.lboiol ). 
the method is likely to be producing distance measure- 
ments for years to come. 
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